o^ 



X 



THE REACTION pp -^ pAK+ NEAR THRESHOLD 

N. Kaiser 

Physik Department T39, Technische Universitat Miinchen, 
D-85747 Garching, Germany 



Abstract 

Q\ I We analyze the recent total cross section data for pp -^ pAK~^ near threshold niea- 

0^ ' sured at COSY. Using an effective range approximation for the on-shell pA S-wave final 

state interaction we extract from these data the combination /C = ^/2\K^f^~+jKt\^ = 
^~r \ 0.38 fm^ of the singlet (Ks) and triplet (Kt) threshold transition amplitudes. We 

Ft , ■ present an exploratory calculation of various (tree-level) vector and pseudoscalar me- 

son exchange diagrams. Pointlike w-exchange alone and the combined {p^,ui, K*~^)- 
exchange can explain the experimental value of /C. The pseudoscalar meson exchanges 
based on a SU(3) chiral Lagrangian turn out to be too large. However, when adding 
^ ■ 7r°-exchange in combination with the resonant vrA^ -^ 511(1650) -^ KK transition and 

^O ! introducing monopole form factors with a cut-off Ac = 1.5 GeV one is again able to 

^T I reproduce the experimental value of /C. More exclusive measurements are necessary 

^\i ■ to reveal the details of the pp — > pAK~^ production mechanism. 

O' 

^ ' PACS: 13.60.Le Meson production - 13.75.Ev Hyperon-nucleon interactions 
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1 Introduction and summary 

With the advent of the proton cooler synchrotron COSY at Jiilich high precision data for 
associated strangeness production in proton-proton collisions, pp — > pAK^, have become 
available in the near threshold region 0, §, §] • These data are of interest in several respects. 
First, they can provide a possibility to test various theoretical models of the strangeness 
dissociation mechanism for the nucleon. Secondly, the cross sections for strangeness produc- 
tion in the elementary pp-collision are an important input into transport model calculations 
of the strangeness production in heavy ion collisions. The latter may provide information 
about the hot and dense state of nuclear matter or the formation of the quark-gluon plasma. 
Various dynamical models have been developed for the reaction pp -^ pAK^ in refs.|^, Q 
which in particular focus on the role of the 5*11 (1650) nucleon resonance decaying into KA. 
In a recent work ^ we have developed a novel approach to pion and eta production in 
proton-proton collisions, pp -^ ppir^ , pmT~^ , pprj , near threshold. In this approach one starts 
from the invariant T-matrix at threshold in the center-of-mass frame which is parametrized 
in terms of one (or two) constant threshold amplitudes. Close to threshold the relative 
momentum of the nucleons in the final state is very small and their empirically known 
strong S-wave interaction plays an essential role in the description of the meson-production 
data. In fact is was found in ref. that in all three cases the energy dependence of the 
total cross section near threshold is completely and accurately determined by the three-body 
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phase space and the on-shell S-wave NN final state interaction (in the case pp — >■ pprj also the 
S-wave ?7A^- interaction was included). Close to threshold the final state interaction can even 
be treated in effective range approximation using the well-known values of the scattering 
lengths and effective range parameters. Note that ref . gives a (partial) derivation of such 
an approach to final state interaction in the context of effective field theory (i.e. using 
only Feynman diagrams). Once one accepts such a phenomenological separation of the 
(on-shell) final state interaction from the full production process, one can extract from 
the total cross section data an experimental value of the constant threshold amplitude 
parametrizing the T-matrix. In the next step a standard Feynman diagram calculation is 
performed for the center-of-mass T-matrix at threshold. It was stressed in ref. |^ that the 
evaluation of the Feynman diagrams has to be done fully relativistically since non-relativistic 
approximations will in general fail to reproduce correctly certain nucleon propagators. The 
source of this problem is the extreme kinematics of the meson-production process with the 



external center-of-mass momentum \p\ ~ J MttIt^^^ being proportional to the square root of 
the nucleon and meson mass. As a major result it was found in ref. that already the well- 
known tree-level (pseudoscalar and vector) meson exchange diagrams lead to predictions for 
the constant threshold amplitudes which agree with the corresponding experimental values 
within a few percent. It also turned out that the short range (p and u) vector meson 
exchange dominates over the long range (vr and r/) pseudoscalar meson exchange for the 
processes pp — > ppn^ , pmT~^ , pprj at threshold. 

The purpose of this work is to present a similar analysis for the kaon production channel 
pp — > pAK^. After defining the threshold T-matrix for pp — ^ pAK^ in terms of a singlet 
(Kg) and a triplet (Kt) transition amplitude and implementing the pA S-wave final state 
interaction in effective range approximation, we extract from the COSY data an experi- 



mental value for the combination /C = J2\Ks\'^ + \Kt\'^ = 0.38 fm^. Next we perform a 
relativistic Feynman diagram calculation of various vector meson (p", u, K*~^) and pseu- 
doscalar meson {n'^,ri,K~^) exchange diagrams with vertices given by a SU(3) symmetric 
(chiral) meson-baryon Lagrangian. We first evaluate these tree diagrams straightforwardly 
from the relativistic SU(3) Lagrangian not introducing ad hoc meson-nucleon form factors. 
The latter are a model-dependent and unobservable concept to account in some sense for 
the finite size of the hadrons involved. It is found that pointlike cj-exchange alone with 
coupling constants given by SU(3) symmetry can well reproduce the experimental value 
of /C = 0.38 fm^. The total [p'^ , uj , K*~^) vector meson exchange leads to /C = 0.45 fm'^ 
which is about 20% too large. Taking into account that possible SU(3) breaking effects 
and the uncertainty of the vector meson baryon coupling constant are of similar size, one 
can argue that the (pointlike) total vector meson exchange is still capable to explain the 
experimental value /C = 0.38 fm^. The pseudoscalar meson {'n''^,ri,K~^) exchange diagrams 
with vertices given by the next-to-leading order SU(3) chiral Lagrangian give rise to rather 
large individual contributions. There is a tendency for cancelation between different types 
of diagrams, but with the SU(3) chiral meson-baryon vertices the effect is not pronounced 
enough. Interestingly, the (pointlike) f^^-exchange alone leads to /C = 0.41 fm^ and if one 
omits the 7r*^-exchange one obtains /C = 0.43 fm^ from the remaining [p^ , u , K*~^ , K~^ , rj)- 
exchange diagrams. It therefore seems that the transition amplitude for irN — > KA as given 
by the SU(3) chiral Lagrangian is not realistic. This is also underlined by the work of ref.[0 
in which these chiral amplitudes have been iterated to infinite orders via a coupled channel 
Lippmann-Schwinger equation and this way a good simultaneous fit of all available low en- 
ergy data for pion (and photon) induced [r], i^)-production could be found. We follow here 
refs.0, Q and add a 7r°-exchange diagram involving the resonant vrA^ -^ S'ii(1650) -^ KA 
transition. We assume the transition strength to be such that this process alone reproduces 



the near threshold data for pp — > pAK^. If we furthermore introduce at each meson-baryon 
vertex a monopole form factor with a cutoff A^ = 1.5 GeV we finally end up with a total 
sum of the S'ii(1650)-excitation graph and the vector and pseudoscalar meson exchange 
diagrams which again is in good agreement with the experimental value /C = 0.38 fm'*. 

It becomes also clear from our analysis that the unpolarized total cross section data for 
pp — > pAK~^ do not provide enough information to distinguish different production mecha- 
nisms. In essence the information given by the total cross section data can be condensed into 



a single number, namely /C = J2\Ks\'^ + \Kt\'^ = 0.38 fm'^. As we have demonstrated here, 
one can find various subprocesses (as e.g. the vector meson exchange, the K~^ -exchange or 
the S'ii(1650)-excitation focussed on in ref.|^) which alone can reproduce this value. More 
exclusive measurements of angular distributions and polarization observables for which one 
does not average out the kinematical complexity of the process pp -^ pAK~^ due to spin and 
the three-particle final state are needed. Such measurements have been started and we 
expect more data to come from COSY in the near future. 

2 Threshold T-Matrix 

The T-matrix for kaon and lamb da-hyp eron production in proton-proton collisions, Pi{p) + 
P2{—p) ^ p + A + K~^, at threshold in the center-of-mass frame reads, 

rr(PP ^ pAir+) = ^ (z ai - ^ ^2 + ^1 X as) ■ J?+ ^ ^(^i + a^) ■ p , (1) 

where p is the proton center-of-mass momentum with \p\ = 861.5 MeV at threshold. The 
spin-operator ai is understood to be sandwiched between the spin-states of the ingoing 
proton pi{p) and the outgoing proton, while (T2 acts between the proton p2{—p) and the 
outgoing lambda. The (complex) amplitude Kg belongs to the singlet transition ^Pq — *• 
^Sqs and the amplitude Kt belongs to the triplet transition ^Pi -^ ^Sis. The factor l/-\/3 
was taken out for convenience since it appears naturally in a calculation employing SU(3) 
symmetry. We follow now the successful approach to pion and eta production of ref . and 
assume the T-matrix to be constant in the near threshold region and the energy dependence 
of the total cross section to be given by the three-body phase space and the (on-shell) pA 
S-wave final state interaction. Since the outgoing proton and lambda have small relative 
momentum one treat their S-wave interaction in effective range approximation, i.e. in terms 
of the singlet and triplet pA scattering lengths and the singlet and triplet pA effective range 
parameters. Experimental evidence [§] and model calculations P suggest that these are 
rather similar for the ^5*0 and ^5*1 pA-states. In this case the unpolarized total cross section 
for pp -^ pAK^ including the pA S-wave final state interaction reads, 

M^MaVs - 4M2 
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^ XiW, M\ Ml) \{W\ ml, s) F^a{W) , (2) 

with e = v^ — M — M\ — ttik the center-of-mass excess energy. M, M\ and itlk denote the 
proton, lambda and (charged) kaon mass. W is the pA invariant mass with values between 
M-|-Ma and the kinematical endpoint y/s—rriK, and A(x, y, z) = x^ +y'^ + z^ ~-2yz —2xz — 2xy 
denotes the Kallen or triangle function. The correction factor from the (equal singlet and 
triplet) pA S-wave final state interaction reads in effective range approximation. 
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We use for the pA scattering length and effective range parameter the values a = 2.0 fm and 
f = 1.0 fm as extracted in ref.[|10| from the Dalitz plot distributions of the pp -^ pAK~^ data 
and the low energy elastic pA scattering cross sections. Using eqs.(2,3) for the total cross 
section, a best fit of the seven COSY data points near threshold p[ leads to the following 
experimental value of the combination 



/C 



2\Ks\^ + \Kt\'^ = 0.38im^ 
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0.68 


1.68 


2.68 


3.68 


4.68 


5.68 


6.68 


^tT [nb] 


2.1 ±0.2 


13.4 ±0.7 


36.6 ±2.6 


63.0 ±3.1 


92.2 ±6.5 


135. ±11. 


164. ± 10. 


<t nb 


2.5 


14.2 


34.0 


60.5 


92.7 


130. 


171. 



Tab.l: Total cross sections forpp 
described in the text. 



pAK~^ . The data are taken from ref.j^ and the fit is 




Fig.l: Total cross sections for pp — > pAK^ as a function of the center- of -mass excess energy 
e = ^/s — M — Ma — rriK- The data are taken from ref.^ and the full line is calculated with 
K, = 0.38 fmf' and pA final state interaction. 

The resulting fit values of o"tot are given in Tab.l and the energy dependent total cross 
section is shown in Fig.l for excess energies e < 7 MeV. Note that the best fit values at the 
two lowest energies e = 0.68 MeV and e = 1.68 MeV lie somewhat outside the experimental 
error band. This may be due to the neglect of the pK~^ Coulomb interaction. A similar 
slight overestimation of the data points closest to threshold was also observed in ref . [Q for 
the reaction pp — > ppn'^. Compared to the processes pp -^ ppn^ , pn7i~^ , pprj studied in ref. [|1 
the pA final state interaction plays here a much less important role. This is also visible from 



Fig.l which shows that the data follow already approximately the pure three-body phase 
space behavior, i.e. o"tot(e) ~ e^- Recently, three additional data points near threshold have 
been measured at COSY 0: atot(e = 8.6 MeV) = (344±41)nb, atot(e = 10.9 MeV) = 
(385 ± 27) nb and crtot(e = 13.2 MeV) = (505 ± 33) nb. These values are well reproduced 
by the present fit which predicts: 261 nb, 383 nb and 518 nb, respectively. We note aside 
that the two data points at higher excess energies (Ttot(e = 55 MeV) = (2.7 ± 0.3) /xb and 
(Ttot(e = 138 MeV) = (12.0 ± 0.4) fjh § come out according to eq.(2) as 3.9 /xb and 11.5 /xb. 
The good agreement for e = 138 MeV should be regarded as accidental, since the same 
experiment has measured a negative A recoil polarization, which can only result from S- 
and P-wave interference terms. 

3 Diagrammatic approach 

In this section we will present an exploratory calculation of various vector and pseudoscalar 
meson exchange diagrams contributing to the threshold amplitudes K^t defined in eq.(l). 
The coupling of vector mesons to baryons is described by the SU(3) symmetric Lagrangian, 

£yB = |l{tr(57^[r^5])+tr(57^5)tr(y^)} . (5) 

The SU(3)-matrices B and V^ collect the octet baryon fields (iV, A, E,S) and the octet 
vector meson fields {p,u,K*,(j)), respectively. The form eq.(5) implies the relations gp^ = 
gv/'2' = QujN /"i and g^N = 0, and we use gv — 6 for the vector meson coupling constant. 
The pseudoscalar meson-baryon interaction is given by the chiral Lagrangians, 

4'i = ^tr(57^[[0, d''<P], B]) + -^tr(5757M{^^</>, B}) + -^tr(5757/.[5'^0, B]) , (6) 
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C% = bDtTiB{x+, B}) + bMB[x+, B]) + botT{BBMx+) , X+ = 2xo - ^{0, {0, Xo}} 

(7) 
with the diagonal matrix xo = diag{rn'^,rn'^,2m1^ — m^). The SU(3)-matrix collects 
the octet pseudoscalar meson fields (vr, K, rf) and f-,, = 92.4 MeV is the weak pion decay 
constant. D ^ 0.75 and F ~ 0.50 are the axial vector coupling constants as determined 
from semi-leptonic hyperon decays. The here most relevant NAK coupling constant is given 
according to eq.(6) as gNAK = — (-D + 3F)(M + MA)/(2v^/7r) = —14.4, which is rather close 
to the empirical value g^AK = —13.2 of ref. |]12|. For the second order chiral Lagrangian C\^ 
only the explicit chiral symmetry breaking terms (linear in the quark masses) are displayed 
in eq.(7). The coefficients 6d,f,o are related to mass sphttings in the baryon octet, the iiN 
(T-term cr7v(0) = (45 ± 8) MeV and the (scalar) strangeness content of the nucleon. We use 
the values bo = 0.066 GeV-\ bp = -0.213 GeV-^ and bo = -0.304 GeV-\ found in ref.[|]. 
At second order there exists in addition a large set of double- derivative terms (~ 9^09y0) 
with a priori unknown coefficients. For the process NN — ^ NNir it was observed in ref.p| 
that the explicit chiral symmetry breaking term is dominant at second order and therefore 
we neglect here the double-derivative terms with unknown coefficients. 

3.1 Vector meson exchange 

The vector meson exchange diagrams contributing to pp -^ pAK^ are shown in Fig. 2. Note 
that according to eq.(5) both the AAp*^ and AS°p° coupling vanish and therefore no p° 
exchange occurs in the right hand diagram where the i^"^-meson is emitted from the proton 
line before the vector meson exchange. 



u},K* 



u,p^ 



Fig. 2: Vector meson exchange diagrams contributing to pp — > pAK^ . 

Straightforward evaluation at threshold of the diagrams shown in Fig. 2 gives for u- 
exchange, 

3g^{D + 3F)mK{4:M - m!) 



H 



K, 



K[^) 



8Mf^{ml + Mm'){2M + m')m' 

?>gl{D + ?,F)mK{m' - M) 
AMf^{ml + Mm'){2M + m')m' 



0.38 fm^ 



-0.07 fm^ 



(8) 



and for p^'-exchange, 



(p) 



K, 



K[^^ 



16M/^(m2 + Mm')(2M + m') 

gl{D + 'iF)mK 
IQMf^iml + Mm'){2M + m') 



m 

am''''' 

2 + -Kr 

. 2 ' 



0.01 fm^ 



m 



8M 



Kp{6 + hKf^ 



-0.12 fm^ . (9) 



Note that we have included the large anomalous tensor-to-vector coupling ratio Kp ~ 6 in 
the ppp°- vertex. For the ppo;- vertex the tensor coupling is known to be small k^^ ~ 0. The 
abbreviation m' stands for m' = M\ — M + mx = 671.0 MeV. Furthermore, one finds from 
i^*"'"-exchange, 



K, 



(K*) 



2Ki'''^ 



SgyFrriK 



2Mf^{rn\, + Mm')m' 



-0.24 fm^ 



(10) 



For the i^*"'"-exchange diagram the intermediate state baryon can be either a A or a S° and 
we neglected the small mass difference M^o — Ma = 77 MeV in eq.(lO). In the absence of any 
empirical evidence, we did not include an anomalous tensor coupling in the pAi^*+-vertex. 
One observes that tu-meson exchange alone, eq.(8), with K,^^^ = 0.39 fm^, well reproduces 
the empirical value /C = 0.38 fm^ extracted from the COSY data. Summing up all vector 
meson (ci;,p°, i^*"'")-exchange contributions one gets /C-^-* = 0.45 fm^ which is about 20% 
too large. Taking into account, that the uncertainty of gy and possible SU(3) breaking 
effects are of similar size one can still argue that (pointlike) vector meson exchange is able 
to explain the experimental value /C = 0.38 fm^. 



3.2 Pseudoscalar meson exchange 

The pseudoscalar meson exchange diagrams contributing to pp —>■ pAK^ are shown in 
Fig. 3. Two types of diagrams are possible, one where the meson rescatters via a chiral 
contact vertex, and another one with a baryon propagating in the intermediate state. 





/ 
/ 








7C^,K+,7] 


/ 
/ 




• 
• 






f 


_ y_^v _ _ 




7t',K+,7] 



Fig. 3: Pseudoscalar meson exchange diagrams contributing to pp —>■ pAK^ 
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The calculation of the rescattering type diagrams in Fig. 3 gives for 7r°-exchange, 
D + F 



(n) 



2ir^) 



16/3 (m2 + Mm') 



A{bD+SbF)iml+ml 



3(m'+m^) =-0.95fm\ (11) 



for fi'^-exchange, 
and for ?7-exchange, 



D + 3F 



if^imj^ + Mm') 



8(&o + bD)ml^ + m' + mx 



-0.87 fm^ 



(12) 



K, 



D-3F 



48/^ 



(m2 + Mm') 



A{bD + 3bF){5m% - 3ml) + 9(m' + m^) 



-0.25 fm^ 



(13) 

From the other diagrams in Fig. 3 with a baryon propagating in the intermediate state one 
finds for vr'^-exchange, 



K, 



M 



2ir^) 



(D + F)mK 



8/3(m2 + Mm'){2M + m') 



4:D{D-F)M+{3D^+2DF+3F'^ 



]m 



for i^^-exchange, 



and for ?7-exchange, 



(K) 



-2Ki^^ 



{D + 3F){D^ + 3F^)mK 



Qfni^K + Mm') 



0.54 fm^ 



(9F2 - D^)mK 



24/3 (m2 + Mm'){2M + m') 



D{AM + m') + 3Fm' 



0.12 fm^ 



0.26 fm^ 

(14) 

(15) 
(16) 



Numerically, the pseudoscalar meson exchange contributions generated by the SU(3) chiral 
Lagrangians eqs.(6,7) are too large, K^ = —1.16 fm^, K^^^^ = —0.25 fm*^. When combined 
with the vector meson exchange terms one would get /C(^+p^) = 1.28 fm^, which is about a 
factor 3 too large. Note that there is some tendency for cancelation between the rescattering 
type diagrams and the other ones, but the effect is not pronounced enough. It is interesting 
to observe that the pointlike f^^-exchange alone eqs. (12,15) gives /C^^-* = 0.41 fm^ and if one 
omits the vr^-exchange eqs. (11,14) one obtains from the remaining contributions /C*-^"'"'^"'"''^ = 
0.43 fm*^. Both these numbers are close to the experimental value JC = 0.38 fm^. It therefore 
seems that the tree level nN -^ KA transition amplitude as given by the next-to-leading 
order chiral Lagrangian eqs. (5,6,7) has some unrealistic features. This is also underlined by 
the recent work of ref . in which these chiral amplitudes have been iterated to infinite order 
via a Lippmann-Schwinger equation and this way a good simultaneous fit of all available 
low energy data for pion (and photon) induced {rj, ii')-production could be found. Instead 
of performing a similar coupled channel calculation we follow here refs.|^, ^ and add a 
further 7r°-exchange diagram involving the resonant vrA^ — » 5*11(1650) —* KA transition. As 
argued in ref.Q] we assume the corresponding transition strength to be so large that this 
process alone can reproduce the near threshold data for pp -^ pAK^ . In this case one 
has the following contribution from S'ii(1650)-excitation to the triplet and singlet threshold 
amplitudes. 



K 



(N*) 



27^1^*) = 0.31 fm^ 



(17) 



We have convinced ourselves that the direct evaluation of the S'ii(1650)-excitation graph 
leads to values of similar size taking into account the empirical ranges of the 5*11(1650)- 
resonance mass and partial decay widths into vrA^ and KA. Let us furthermore introduce a 



monopole form factor at each (off-shell) meson-baryon vertex with a cut-off Ac = 1.5 GeV. 



Such values of the cut-off Ac are typically used in one-boson exchange models [|T3| of the 
NN-interaction for both the pion- and vector-meson-nucleon vertex. This modification of 
the diagrams brings a reduction factor (1 + Mm'/A^)~^ for the vector and pseudoscalar 
meson exchange graphs. Note that in the case of the S'ii(1650)-resonance contribution 
the form factor effect is already included in the number given in eq.(17). Summing up 
all the contributions due to S'ii(1650))-excitation as well pseudoscalar and vector meson 
exchange (including the form factor) one finds KI° = —0.31 fm^ and K^^°^^ = —0.18 fm^. 
The resulting value /C(*°*) = 0.40 fm^ is again in good agreement with the empirical value 
/C = 0.38 fm^ 

Evidently, the main lesson to be learned from the present exploratory calculation is 
that there are in fact various subprocesses (like pointlike vector meson or i^+-exchange or 
the S'ii(1650)-resonance excitation) which alone can explain the near threshold data for 
pp — > pAK~^. With reasonable assumptions on the various coupling strengths and the cut- 
off Ac entering the (unobservable) meson-nucleon form factor one finds that also the total 
sum of many processes is able to reproduce the near threshold data for pp -^ pAK~^. This is 
possible because of cancelations between terms of different sign and because of the freedom 
the to shift strength between the singlet (Kg) and triplet (Kt) threshold amplitude. Only 
more exclusive data (like angular distributions and polarization observables) can help to 
distinguish different pp -^ pAK~^ production mechanisms. 

Finally, we like to comment on the recently measured process pp — > pTi^K'^. It was found 
experimentally |]TT| that the corresponding total cross sections near threshold are about a 



factor 30 smaller than those for pp — > pAK^ (at equivalent excess energies), cu- meson 
exchange or i^+-exchange might offer an explanation for this suppression via the small 
SU(3)-ratio Qnhk/ Qnak = — V^/9 ^ —0.19. An explicit calculation of the a;-exchange and 
if "'"-exchange (using now m' = M^o — M + rriK = 748 MeV in eqs. (8,12,15)) gives for the 
pTi^ K~^ -channel values of JC^J = 0.06 fm^ and /C^o ~ ^-^^ im^, i.e. a suppression factor of 
about 40 or 20 for the total cross section. However, such considerations may be too simplistic 
in the light of possible strong coupled channel effects 0] between KA- and KS-states. 
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